Stable suspension and dispersion-induced transitions 
from repulsive Casimir forces between fluid-separated cylinders 



oo 

o 

O 
> 

o 
oo 



43 

Oh 
i 

cr 



> 



O 
OO 

O 
> 

X 



Alejandro W. Rodriguez, 1 J. N. Munday, 2 J. D. Joannopoulos, 1 Federico Capasso, 3 Diego A. R. Dalvit, 4 and Steven G. Johnson 5 

Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

Department of Physics, Harvard University, Cambridge, MA 02139 

School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02139 

4 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

5 Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139 

We numerically demonstrate a stable mechanical suspension of a silica cylinder within a metallic cylinder 
separated by ethanol, via a repulsive Casimir force between the silica and the metal. We investigate cylinders 
with both circular and square cross sections, and show that the latter exhibit a stable orientation as well as a 
stable position, via a method to compute Casimir torques for finite objects. Furthermore, the stable orientation 
of the square cylinder undergoes a 45° transition as the separation length-scale is varied, which is explained as 
a consequence of material dispersion. 



The Casimir force, arising from quantum fluctuations of 
the electromagnetic field, was first described as an attractive, 
monotonically decaying force between metallic plates fT), 
but repulsive interactions can arise in special circumstances, 
e.g. involving fluid-separated asymmetric plates [2|. It has 
been proposed that repulsive Casimir forces between fluid- 
separated objects can lead to stable mechanical equilibria, and 
hence frictionless static bearings or other interesting passive- 
suspension devices Q]|4]. However, previous calculations and 
experiments involving fluid-separated objects have been re- 
stricted to geometries involving parallel plates or approxima- 
tions thereof |5| (similar to work on air-separated metals, as 
reviewed in several recent papers, e.g. Ref. 6). Here, using 
recently developed numerical techniques [7], we present the- 
oretical calculations of Casimir forces/torques between fluid- 
separated objects with finite square and circular cross sec- 
tions (Fig. [T] insets) that rigorously demonstrate stable po- 
sitional/orientational Casimir equilibria. (Vacuum-separated 
perfect-metal circular cylinders were shown to exhibit an un- 
stable positional equilibrium [8 1.) In the case of square cross 
sections, a surprising result is obtained for the rotational equi- 
libria: the stable orientation changes by 45° depending on the 
lengthscale, which we show to be a consequence of material 
dispersion. In particular, for certain fluid-separated materi- 
als, the frequency dependence of the permittivity e (material 
dispersion) causes the Casimir force to switch from repulsive 
to attractive at some critical separation [2J; this leads to the 
orientation transition described here, and may produce other 
lengthscale-based qualitative transitions in future geometries. 
We present a method to compute Casimir torques on finite ob- 
jects (in contrast to previous uncontrolled approximations [9| 
or analytical results for planar cases (3] [TO] [TTJ [TJl [13] [141 ): 
we supplement the accurate results with a heuristic model 
based on the proximity force approximation (PFA) that turns 
out to capture qualitative behaviors of the orientation transi- 
tion, and therefore provides some simple insight. 
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A repulsive Casimir force arises between parallel plates of 
£1.2 separated by a fluid of ey, if £i(?£) < £/(*£) < £2(1%) 
for a sufficiently wide range of imaginary frequencies £ = 
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FIG. 1: Casimir force F x in the x-direction, per unit z -length, on a 
Si02 cylinder suspended within a perfectly-metallic cylinder (inset), 
separated by fluid (ethanol), as a function of the ^-displacement from 
equilibrium (eccentricity) d [in units of a — 0.5(D — s)] for both 
circular (solid-line) and square (dashed-line) cylinders, d — is seen 
to be a stable equilibrium. 



Imw (3). Three such materials are silica (SiC^) and metal 
(Au) separated by ethanol [4|, as discussed below. We stud- 
ied the three-dimensional constant cross-section geometries 
shown in the inset of Fig. [T[ square or circular SiC>2 cylin- 
ders of diameter s surrounded by a metal cylinder of diameter 
D, separated by ethanol. For computational ease, we use a 
perfect metal [e(i£) — ► 00] for the outer cylinder, but we use 
experimental w-dependent e for the SiC<2 and ethanol. In such 
cases, with object sizes and curvatures comparable to their 
separation, approximations of the Casimir force as a pairwise 
attraction between surfaces are not valid and their qualitative 
predictions (e.g. stability) may be incorrect 171 fT31 [161 . 

Figure [T] shows the force per unit z-length on the inner 
cylinder as a function of the displacement d from equilibrium 
in units of a, where a = 0.5(1? — s) = 0.0955 /im so that 
d/a = ±1 for touching surfaces, for parameters s/D = 0.25. 
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FIG. 2: Casimir torque t z , per unit z-length, on the inner square rod 
of eccentric square cylinders as a function of the angle 9 with respect 
to the :r-axis (see right insets) for two material choices as shown in 
the top insets. Error bars are estimates of the effect of finite grid 
resolution. 



It demonstrates a stable equilibrium for both square (solid 
black) and circular (dashed grey) cross sections. The com- 
putational method is based on integration of the mean elec- 
tromagnetic stress tensor (valid even for fluids 03)) evalu- 
ated in terms of the imaginary-frequency Green's function via 
the fluctuation-dissipation theorem [7|. Additional details are 
provided below. 

Having demonstrated positional stability (Fig. fTJ, we now 
explore an unusual effect arising from material dispersion. 
The square-cylinder geometry (Fig. [2] inset) must, by sym- 
metry, exhibit two equilibrium orientations (9 = 0° and 
9 = 45°). A priori, it is not clear which of the two configu- 
rations is stable. To determine this, Fig. [2] plots the Casimir 
torque per unit z-length on the inner square, as a function 
of the rotation angle 9 between the inner and outer squares. 
In addition to computing the torque for the SiC»2-ethanol- 
perfect-metal case, we also analyze the torque for vacuum- 
separated perfect-metal cylinders, in which case the forces 
are purely attractive and there is no stable positional equilib- 
rium |8|. 

The resulting torque per unit z-length is shown in Fig. [2] 
for both the repulsive (solid, filled-square-lines) and attrac- 
tive (solid, open-square-line) cases, as 9 is varied from —45° 
to 45°. The stable orientation for both cases is 9 = 45°, a 
surprising result considering the difference in the sign of the 
force. A stable orientation of 45° in the perfect-metal case is 
not so surprising, since an attractive force that decreases with 
separation should intuitively (in the heuristic picture of pair- 
wise attractions between surfaces) favor an orientation where 
the surfaces are as close as possible, pulling the corners of the 
inner square towards the outer surfaces. On the other hand, 
a repulsive force should intuitively push the surfaces as far 
apart as possible, which would suggest a 0° stable orienta- 
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FIG. 3: Derivative of the Casimir torque in the ^-direction with re- 
spect to 8, d,T z /d6, in units of he/ a 2 , evaluated at 8 ~ 0°, as a 
function of lengthscale a = 0.5(D — s). The solid line is a fit to 
the exact Casimir torque, computed by our numerical method (solid 
squares), and the dashed line is the torque as computed by a PFA 
approximation. Both curves display a change (dispersion-induced 
transition) in the stable orientation of the square. 



tion. The reason for this apparent contradiction between in- 
tuition and the numerical results lies in the effects of mate- 
rial dispersion — as explained below, esi ~ £cth (and hence the 
force) switches sign at some lengthscales (some I/O- F° r 
example, if we neglect the material dispersion of SiC<2 and 
ethanol, and use only the £ — » dielectric constants, the 
torque as a function of 9 (shown as the dashed black line in 
Fig. Ell indeed exhibits a stable orientation at 0°, as expected 
from the pairwise intuition: in this case, the force is repulsive 
at all lengthscales. 

Since 9 = 0° is unstable for the a = 0.0955 /im fluid 
case and stable for a — > oo (equivalent to £ — > 0), there 
must be a transition at some critical intermediate lengthscale 
a c . One way to determine a c is to calculate the derivative of 
the torque at 9 = 0° as a function of a (noting that the di- 
mensionless torque to 2 /he and its derivative go to a nonzero 
value as a — ► oo), and to look for a change in the sign. This 
derivative is plotted as the solid line in Fig. [3] and displays 
the expected transition from unstable (d,T z /d9 > 0, left) to 
stable (d,T z /d9 < 0, right) at a c fw 0.1 /im, a consequence 
of material dispersion. A better understanding of this tran- 
sition can be gained by inspecting a simple heuristic (PFA). 
PFA is only an ad hoc model, in which the force on each 
point of the surface (and hence the torque) is treated as sim- 
ply the parallel-plate (Lifshitz) force between fluid-separated 
half-spaces (with "lines of interaction" perpendicular to the 
inner square, as depicted in the left inset of Pig . [31 . This sim- 
ple model turns out to capture some qualitative features of the 
orientation transition, as shown by the dashed line in Fig. [3] 
although it is of course quantitatively incorrect iflTTl . 

PFA provides an explanation for why orientational stability 



need not coincide with positional stability. Stability, in gen- 
eral, arises from the competing interactions of the inner and 
outer surfaces, e.g. in the 45° orientation for repulsive inter- 
actions, the nearest surface "pushes" the corner away, while 
the other surface "pushes" the comer back. Which of these 
competing forces dominates depends on their sign and power 
law, but in the torque r = r x F case, there is an additional 
effect (from the r x dependence): the force "from" the near- 
est surface elements is parallel to the radial direction and does 
not produce a torque, whereas it does contribute to the total 
force. Therefore, the net force and torque can respond differ- 
ently to the same material dispersion, because the lengthscales 
of the dominant contributing separations differ. More gener- 
ally, in the exact model where the force cannot be decomposed 
into additive contributions, the dominant imaginary-frequency 
contribution can differ between torque and force. 

There is another interesting feature in Fig. B] for a > a c , 
the derivative of the torque is nonmonotonic, decreasing in 
magnitude towards a — ► oo. This may seem counter-intuitive 
because the e contrast is maximum for £ — ► (Fig. HJ inset), 
which would seem to predict greater forces. However, be- 
cause the net torque arises from a competition between nearby 
and faraway surfaces, the magnitude depends not only on the 
force but also on the power law: if the force decreases more 
rapidly with distance, then faraway surfaces contribute less 
and the net torque is larger. This is precisely what is happen- 
ing here: at intermediate a, material dispersion acts to increase 
the power law [2 1 compared to the a — > oo case where disper- 
sion plays no role. Preliminary work shows similar effects in 
the positional stability. This also explains the relative torque 
magnitudes for a — ► oo and a = 0.0955 fjxn in Fig. [2] 

Since the qualitative behavior of the orientation transition 
is captured by PFA, it makes sense to examine a simple one- 
dimensional model in detail, depicted on the inset of Fig. [4] 
a SiC>2 slab of thickness s, separated from two metallic half- 
spaces, a (surface-to-surface) distance a = 0.5(1? — s) from 
the SiC>2 surface, with ethanol (fluid) in between. This ge- 
ometry is analyzed using a generalization of the Lifshitz for- 
mula |18|. Figure HI shows the force per unit area on the in- 
ner slab (again using s/D = 0.25), as a function of the dis- 
placement d, for two lengthscales: a = 0.0955 /im (red), and 
a = 0.6 /im (black). There is a stable equilibrium at d = 0, 
just as for the two-dimensional case in Fig. W\ In the case 
of perfect-metal half-spaces, the force changes sign at a crit- 
ical distance d c that depends on the lengthscale a. However, 
if we use a realistic metal (Au), shown by the dashed line 
for a = 0.0955 /im, this transition disappears and the force 
is always repulsive. These features are explained by the in- 
set of Fig. El which shows e(i£) for SiC>2, ethanol, and Au. 
The key point is that the contributions to the Casimir force 
come primarily from imaginary "wavelengths" 2n/£ larger 
than some lengthscale set by the separation, while very short 
wavelengths (large £) on this scale are exponentially cut off 
in the force integral [19]. Thus, for large lengthscales a, the 
force is dominated by small-£ contributions where e s iii ca < 
Methanol < £metah giving a repulsive force and a stable equi- 
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FIG. 4: Casimir force F x in the a>direction, per unit area, be- 
tween a planar S1O2 slab suspended between two perfect-metal plates 
(solid red and black lines) or gold half-spaces (dashed red line), 
separated by a fluid (ethanol), as a function of the dimensionless 
x-displacement d from equilibrium. The force is plotted also at 
two lengthscales a = 0.5(D — s) = 0.0955 /im (red lines) and 
a — 0.6 /im (black line). The solid blue line shows the same quan- 
tity for the case of a lithium-niobate slab and a gold half-space, for 
a = 0.5 /im. (Insets:) Top inset: e(i£) for Si02 (black), LiNbOa 
(solid and dashed blue), ethanol (orange) and gold (red) as a func- 
tion of imaginary frequency £. Bottom inset: schematic of geometry. 



librium. On the other hand, for small lengthscales and separa- 
tions, large-£ contributions become more important, for which 
Esiiica > Ecthanoi (for £ > 2.3 c//itn, marked by the verti- 
cal dashed line in the inset of Fig. HI), leading to attractive 
forces. For the case of perfect-metal cylinders, these attrac- 
tive contributions are large enough to flip the sign of the to- 
tal force for small separations, whereas for Au the attractive 
high-£ contributions are suppressed by the diminishing e(i£) 
of the Au and the force remains repulsive. While the transition 
is absent in the SiC>2 -ethanol- Au configuration, we have cal- 
culated the force between ethanol-separated lithium-niobate 
LiNbOa and Au plates and found a sign transition at separa- 
tions ~ 0.05 /im. This is illustrated in Fig.HJ(solid blue line); 
because LiNbCVs anisotropy greatly complicates the mod- 
eling, we used a PFA approximation based on summing the 
Lifshitz forces between semi-infinite Au and LiNbO,3 slabs 
(neglecting the finite LiNbO.3 thickness). The two principal 
values of the e tensor of LiNb03 [20] are plotted in the in- 
set. There are also transitions for ethanol-separated barium- 
titanate and calcite plates at separations ~ 0.01 /im [3]. Since 
these two cases are much more difficult to compute (the outer 
cylinder is not perfect metal and LiNbOs is anisotropic), we 
focused on the Si02-ethanol-metal case here, which has sim- 
ilar qualitative behaviors. 

The experimental es of SiC>2 and ethanol were fit to 
a standard multiple-oscillator model [21 1 accurate from 
infrared to ultraviolet wavelengths: e(i£) = 1 + 
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1221 : efhanol(7V = 2) u. 



The parameters we used are | 

{6.6, 114} x 10 14 Hz and C n = 
{23.84, 0.852}; Si0 2 (JV = 3) uj n = {0.867, 1.508, 203.4} x 
10 14 Hz and C n = {0.829,0.095,1.098}. We model the 
dielectric constant of Au by a Drude model e(i£) = 1 + 
Wp/£(£ + 7), where u p = 1.367 x 10 14 Hz and 7 = 
5.320 x 10 13 Hz |20). The Casimir force F is computed as an 
integral F ~ L d£,$(T)dA [71, where the mean stress ten- 
sor (T) is computed at each position and frequency from the 
Green's function as described in Ref. 7. The Casimir torque 
r, was computed in the same framework, with a minor mod- 
ification, t ~ J Q d£<jij)r x ((T)dA), also proposed in our 
earlier work Q. 

We are hopeful that these phenomena will be amenable to 
experiment. There are some experimental advantages over 
planar geometries: the concentric configuration is a stable 
equilibrium, static charges on the outer metallic cylinder are 
screened in the interior, and the equilibrium is easily distin- 
guished from electrostatic effects (which cannot produce sta- 
bility [ 15 1). (If we inverted the geometry, to have the Si02 on 
the outside and the metal on the inside, then the larger static e 
of the fluid compared to the SiC<2 would lead to a classical sta- 
ble equilibrium for a charged inner cylinder Ifl"5ll23ll24l .') We 
expect that similar stability will be obtained for real metals, 
as in our Au calculation in Fig. HI and there are several mate- 
rials that exhibit repulsive-attractive transitions (unlike Si02- 
ethanol-Au) that should display the orientation transition. We 
have also uncovered an intriguing question to explore in future 
work: in the PFA heuristic, it appears that the critical length- 
scale at which the orientation transition occurs is smaller than 
the corresponding positional transition lengthscale (changing 
from stable to unstable suspension), whereas the exact calcu- 
lations give the opposite result. 
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